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SUMMARY 

Membrane analysis of very thin pressurized spheroid 
shells composed of flat gores led to a set of nonlinear dif- 
ferential equations for three membrane displacements. An 
approximate solution by the perturbation technique is es- 
tablished. As an example, the approximate shape of Echo II 
is computed. 
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MEMBRANE ANALYSIS OF PRESSURIZED 
THIN SPHEROID SHELLS 
COMPOSED OF FLAT GORES, 

AND ITS APPLICATION TO ECHO II 

by 

Hossein Bahiman 
Goddard Space Flight Center 


INTRODUCTION 

To compute the shape of a pressurized thin spheroid shell composed of many identical flat 
gores requires solution of five nonlinear partial differential equations of equilibrium and boundary 
conditions for three displacements of a gore. The position of each point on the gore could be de- 
termined as a function of two independent variables, namely the point T s distance from the center 
line of the gore and latitude from the equator of the spheroid. For an approximate solution to these 
nonlinear partial differential equations each displacement is expanded into a power series in the 
point's distance from the gore center line, such that the series coefficients are at most functions 
of the latitude. By equating to zero the terms of lowest power in x in each equation, five new non- 
linear differential equations for the six coefficients (two for each displacement) are obtained. The 
approximate solution of these equations by the perturbation technique led to three displacement 
functions of the gore, and one unknown constant which was evaluated by applying the variational 
method to the total elastic strain and potential energy of the system. 

In this analysis, it is assumed that the displacement of the spheroid due to gravity is negligible 
compared to that caused by the constant internal pressure. It is also assumed that the material of 
the spheroid is homogeneous and isotropic and it obeys Hooke's law. 

Finally, this method of analysis is applied to Echo n, the 135-foot diameter passive communi- 
cation satellite. However, since the 3-ply aluminum -Mylar -aluminum material of the satellite 
does not satisfy the aforementioned assumptions, the analytical results obtained here are expected 
to be only an approximation of the actual case. 

THEORY 

To determine the shape of a thin spheroid shell composed of identical flat gores it is sufficient 
to compute the true shape of a gore in the pressurized state. To do this, it is assumed that the 
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spheroid material is homogeneous, isotropic, and linearly elastic. It is also assumed that the 
effect of gravity on the spheroid shape is negligible compared with that of the constant internal 
pressure. 

Let rectangular Cartesian coordinate axes (Figure l)be drawn so that the origin 0 is located 
at the center of the balloon, the y axis is along the polar axis, and the z axis is extended through 
the center of a gore. Draw x' y' z' axes from a point P of the gore in such a way that the x' axis 
will be parallel to the x axis, the z ' axis will be in the direction OP, and the y' (or 0 ) axis will be 
tangent to the meridian at point P. 


In the pressurized state of equilibrium, the resultant of all the forces acting on an element of 
the shell in the direction x', 6 and z' vanishes. Let N x , N 09 and N x<9 be normal and shear tractions 
per unit length of membrane, 9 be the angle between the position vector and the equatorial plane, 
and p be the internal pressure of the spheroid (Figure 2). Let u, v, and w be displacement com- 
ponents at a point P in the directions x', y', and z' respectively. Then (Reference 1) 
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Figure 1 —Coordinate axes and displacement components. 


Figure 2— Forces acting on an element of a gore. 
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The elastic strains are given by: 
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where K = Et/(l - ix 2 ) is the membrane stiffness, and E, t, and \x are Young’s modulus, the shell 
thickness, and Poisson’s ratio, respectively. Substitution of Equations 2 into Equations 3 leads to 
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The boundary conditions are: 


N x £ - 0 On x - ra cos 9 I 

V ( 6 ) 

u - aw cos 9 - av sin 9 On x = ra cos 9 \ 

along the seam, where a = 77 /n. The first expression indicates that because of symmetry shear 
stress vanishes along the seam where any two adjacent gores meet; and the second expression 
implies that the points along each seam move in the meridional plane. Since the coefficients of 
only the first two terms of the power-series expansion of displacements will be determined by this 
method, continuity of slope in the plane normal to the seam cannot be used as a third boundary 
condition. The use of such a boundary condition would lead to an erroneous solution, namely, that 
the cross section of the spheroid with any plane normal to its polar axis is a true circle, no matter 
what the magnitude of the internal pressure. Instead, a virtual displacement is given to the pres- 
surized spheroid in the state of equilibrium and change of the total elastic strain and potential 
energy of the spheroid and the load is equated to zero. 

Let r be the radius of the centerline of a gore before any load is applied, and a, b, c, d, f and 
g be functions of the variable angle 9 . Displacements w, u, and v may be expanded into power series 
as follows: 
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Substitution of Equations 7 and the last of Equations 4 into Equations 5 and 6 leads to 


(?) 


6d + 



dg 




da 

db 1 

- jJL d 2 c 


d 2 a 

+ 4b 2 


(i 

f M) 

d e + 

2/jb + 

(1+M) 

d9 

d# + “ 


<1-M) 

b d& 2 

= 0 

1 +/X 

dc 






da 

dii 

da da 

d 2 a 

db 



2 

dB 

+ ( 

1 

g + (1 

-3/-0 

b 

d# 

+ d# 2 

+ A6 + d6 

d# 2 

2/^a 

= 0 , 


fdf 


1 / 

f da \ 2 


- 

( 


d 2 a \ 



1 

daN 2 

pr 

[d 9 ' 

f a + 

2 ' 

<d9) 

+ fJUC “ 

2/J.ab 

V 

L - 

d0 2 t 

2bc - 2/xb 

- 2/^ab 

~ Mb ( 

,d e) = 

IT 




> 


(8) 


dc da 

d0 + 2g + 2b de - 0 , 


c - a + ft an# 


(b - d) a 2 cos 2 6 


J 


4 



I 


Neglecting small terms of higher order in Equations 8 leads to 
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Since n, the number of gores, is assumed to be large, a 2 is very small. Thus the perturbation 
technique can be applied to the functions a, b, c, d, f, and g in the following fashion: 
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Substitution of Equations 10 into Equations 9 yields a set of five nonlinear differential equations in 
terms of a 0 , b 0 , c 0 , d 0 , f 0 , g 0 , a 1? b 1? c 1? d 1? etc. containing even powers of a. Terms containing 
each power of a in each differential equation must independently vanish. The terms independent 
of a must vanish in the following way: 
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The solution of the above equations could be given as the following: 
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The coefficients of a 2 in Equations 9 must vanish in the following fashion: 
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The solution of Equations 13 could be written as: 
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in which the constant coefficients are given by 
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Knowing the value of Poisson’s ratio, one can easily evaluate in an orderly manner y , f 1? A 1? c 1? 
a 2 , C 2 , D t , and g x as a function of the Constant B. 

To evaluate constant B, the method of virtual displacement may be utilized. The elastic strain 
energy of one quarter of a gore is given by 
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Substitution of stress and strain values in Equation 16 yields 
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In a virtual displacement corresponding to change of b by SB the work done by the internal pressure 
for a quarter of a gore is given by 
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in which w is the total work done by the internal pressure for a quarter gore, and w is the radial 
displacement given by Equations 7. Hence, the Constant B could be determined from 
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Thus for a given pressure p, Equation 19 in conjunction with Equations 17, 7, 10, 12, 14 and 15 
furnishes B, which upon substitution into Equations 15, gives coefficients A 1? A 2 , etc. Substitution 
of these coefficients into Equations 14 in conjunction with Equations 7, 10 and 12 gives the three 
displacement functions u, v, and w. 


CALCULATION FOR ECHO II 


The 135 -foot diameter Echo II is composed of 106 identical gores with a maximum width of 
4 feet which occurs in the equatorial plane of the balloon. The 3 -ply material of the balloon is 
composed of 0.00035 -inch thick Mylar sandwiched between two layers of aluminum each 0.00018- 
inch thick. The adjoining gores are located edge-to-edge with one-inch wide tape of the same 
material sealing the seams lengthwise. When separation of the two half-canisters takes place in 
the initial phase of orbiting, the residual air and water vapor inside the folded balloon— which is of 
the order of 1 mm Hg—inflates the satellite. While initially flat, each gore forms a narrow trans- 
verse portion of a half circular cylindrical surface 135 feet in diameter. Thus the initial shape of 
the balloon is a spheroid with 0.356-inch maximum deviation from a sphere which occurs in the 
equatorial plane. 

Because of solar radiation the skin temperature increases and as a result inflation material 
sublimes inside the open balloon, building up internal pressure until a maximum of 225 ^ Hg is 
reached; this maximum pressure corresponds to a membrane stress of 1.7 lb/in. 

Although the material of the balloon is rather nonhomogeneous and anisotropic and does not 
quite obey Hooke's law, for the purpose of approximation it may be assumed that the requirements 
of homogeneity, isotropy and linear elasticity necessitated by this analysis are met. The average 
value of K corresponding to a skin traction of 1.70 lb /in. (given by Tests 1, 3, 4, and 6 of Refer- 
ence 2) is 1050 lb/in. 
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Since nominal Young 1 s modulus for aluminum is 10. 5 X 10 6 psi and for Mylar is 0.55 X 10 6 psi, 
therefore it is expected that aluminum should carry most of the load. The nominal value of k for 
aluminum layers is given by 

Et 

K_ ~ T - 4240 lb/ inch , 

l - fji 2 

in which m = 0.33 has been inserted. This value of K n is much larger than the experimental value 
K = 1050 lb/in. Therefore, the actual Young's modulus and perhaps Poisson's ratio for thin layers 
of aluminum are smaller than the nominal ones. 

Following the steps suggested in the last paragraph of the previous section, one could get an 
equation for B whose coefficients are functions of parameter m = pr/Ka 2 . For different values of 
maximum pressure B is computed and listed in Table 1. Hence forp = 225^ Hg (or 0.000435 lb/ in 2 ) 

0.000435 x 810 

m = ~ — = 3.81 , 

1050 (ijjg) 


Table 1 

Values of B versus m. 
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2.0 

-.294 

4.0 

-.320 

6.0 
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8.0 
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.1 

-.043 

2.1 

-.298 

4.1 

-.320 

6.1 

-.320 

8.1 

-.317 

.2 

-.079 

2.2 

-.300 

4.2 

-.320 

6.2 

-.320 

8.2 

-.317 

.3 

-.109 

2.3 

-.303 

4.3 

-.320 

6.3 

-.319 

8.3 

-.317 

.4 

-.135 

2.4 

-.305 

4.4 

-.320 

6.4 

-.319 

8.4 

-.317 

.5 

-.157 

2.5 

-.307 

4.5 

-.320 

6.5 

-.319 

8.5 

-.317 

.6 

-.176 

2.6 

-.309 

4.6 

-.320 

6.6 

-.319 

8.6 

-.316 

.7 
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2.7 
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4.7 
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6.7 
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8.7 

-.316 

.8 

-.207 

2.8 

-.312 

4.8 

-.321 

6.8 

-.319 

8.8 

-.316 

.9 

-.220 

2.9 

-.313 

4.9 

-.321 

6.9 

-.319 

8.9 

-.316 

1.0 

-.231 

3.0 

-.314 

5.0 

-.321 

7.0 

-.319 

9.0 

-.316 

1.1 

-.241 

3.1 

-.315 

5.1 

-.321 

7.1 

-.318 

9.1 

-.316 

1.2 

-.250 

3.2 

-.316 

5.2 

-.321 

7.2 

-.318 

9.2 

-.316 

1.3 

-.258 

3.3 

-.317 

5.3 

-.321 

7.3 

-.318 

9.3 

-.316 

1.4 

-.265 

3.4 

-.318 

5.4 

-.320 

7.4 

-.318 

9.4 

-.316 

1.5 

-.272 

3.5 

-.318 

5.5 

-.320 

7.5 

-.318 

9.5 

-.315 

1.6 

-.277 

3.6 

-.319 

5.6 

-.320 

7.6 

-.318 

9.6 

-.315 

1.7 

-.282 

3.7 

-.319 

5.7 

-.320 

7.7 

-.318 

9.7 

-.315 

1.8 

-.287 

3.8 

-.319 

5.8 

-.320 

7.8 

-.317 

9.8 

-.315 

1.9 

-.291 

3.9 

-.320 

5.9 

-.320 

7.9 

-.317 

9.9 

-.315 
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which corresponds to B = - 0.319 in Table 1. Thus, the maximum differential radial displacement 
at the midgore and at the seam which occurs in the equatorial plane is given by 
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x = 0 
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ii 
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0.227 in 


This corresponds to a maximum radial difference 


A = 0. 356 - 0. 227 = 0. 129 in 

in which 0.356 in * is in its initial unpressurized value. 


CONCLUSIONS 

The maximum differential radial displacement of a pressurized spheroid composed of identical 
flat gores, in a zero gravity field, as computed by this technique is believed to be a very good ap- 
proximation of the true one, provided that the material is homogeneous, isotropic, and linearly 
elastic. 


(Manuscript received February 3, 1965) 
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